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Abstract 

We study Stirling permutations defined by Gessel and Stanley in 
We prove that their generating function according to the number 
of descents has real roots only. We use that fact to prove that the 
distribution of these descents, and other, equidistributed statistics on 
these objects converge to a normal distribution. 

1 Introduction 

In [6] Ira Gessel and Richard Stanley defined an interesting class of multiset 
permutations called Stirling Permutations. Let Qk denote the set of all 
permutations of the multiset {1, 1, 2, 2, • ■ ■ , n, n} in which for all i, all entries 
between the two occurrences of i are larger than i. For instance, Q2 has three 
elements, namely 1122, 1221, and 2211. It is not difficult to see that Qn has 

1-3 (2n — 1) = (2n — 1)!! elements. Gessel and Stanley then proved many 

enumerative results for these permutations and showed several connections 
between these and other combinatorial objects, such as set partitions. 

Counting Stirling permutations by descents, the authors of [6j found a 
recurrence relation similar to the recurrence relation known for classic per- 
mutations. In this paper, we will continue in that direction. First, we show 
the simple but interesting fact that on Q„ the descent and the plateau statis- 
tics, to be defined in the next section, are equidistributed. Then we prove 
that for any fixed n, the generating polynomial of all Stirling permutations 
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in Qn with respect to the descent statistic has real roots only. This is analo- 
gous to the well-known case (see Theorem 1.33 of [Ij) of classic permutations, 
namely the result that all the roots of Eulerian polynomials are real. Finally, 
we apply a classic result of Bender to use this real roots property to prove 
that the descents of Stirling permutations in are normally distributed. 

2 Stirling Permutations and Real Zeros 

Let q = 0102 • • • 02n G Qn be a Stirling permutation. Let the index i be 
called an ascent of g' if i = or Oj < Oj+i, let i be called a descent of q 
if i = 2n or Oj > Oj+i, and let i be called a plateau of g if Oj = Oj+i. It 
is obvious that the ascent and descent statistics are equidistributed, since 
reversing an element of Qn turns ascents into descents and vice versa. It is 
somewhat less obvious that the plateau statistic is also equidistributed with 
the previous two. This fact, and a reason for it, are the content of the next 
proposition. Note that its first identity, ([T]), was proved in [6]. 

Proposition 1 Let Cn,i be the number of elements of Qn with i descents. 
Then for all positive integers n,i > 2, we have 

Cn,i = iCn-i,i + (2n - i)Cn-l,i~l. (1) 

Similarly, let Cn,i be number of elements of Qn with i plateaux. Then for all 
positive integers n,i >2, we have 

Cn,i = iCn-l,i + {'^n-i)Cn-l,i-l. (2) 

In particular, since Ci^i = ci^i = 1 and Ci^o = ci,o = 0, the identity 

Cn,i — Cn^i (3) 

holds. 

Proof: There are two ways to obtain an element of Qn from an element 
p G Qn-i by inserting two copies of n into consecutive positions. Either p 
must have i descents, and then we insert the two copies of n into a descent, 
or p has i — 1 descents, and then we insert the two consecutive copies of n 
into one of the (2n — 1) — (i — 1) = 2n — i positions that are not descents. 
The argument proving ([2]) is analogous. O 
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Corollary 1 On average, elements of Qn have (2n+l)/3 ascents, (2n+l)/3 
descents, and (2n + l)/3 plateaux. 

Proposition [1] enables us to prove a strong result on the roots of the 
polynomials Cn.i-c^- The method we use follows an idea of H. Wilf 

(0, [I] Theorem 1.33) who used it on classic permutations. 

Theorem 1 Let Cn{x) = Yl7=i ^n,ix'^ ■ Then for all positive integers n, the 
roots of the polynomial Cn{x) are all real, distinct, and non-positive. 

Proof: For n = 1, one sees that Ci(x) = x, and the statement holds. For 
n = 2, one sees that 6*2(2;) = 2x^ + x = x{2x + 1), and so the statement 
again holds. 

For n > 3, recurrence relation ([T|) implies 

Cn{x) = {x- x2)C;_i(x) + (2n - l)xCn-i{x) (4) 

as can be seen by equating coefficients of x*. The right-hand side is similar 
to the derivative of a product, which suggests the following rearrangement 

Cn{x) = X(l - X)2" A ((1 _ x)l-2-C„_i(x)) . (5) 

Let us now assume inductively that the roots of C„_i(x) are real, distinct 
and non-positive. Clearly, C„(x)) vanishes at x = 0. Furthermore, by Rolle's 
theorem, ([5]) shows that C„(x) has a root between any pair of consecutive 
roots of Cn-i{x). This counts for n — 1 roots of Cn{x). So the last root 
must also be real, since complex roots of polynomials with real coefficients 
must come in conjugate pairs. 

There remains to show that the last root of C„(x) must be on the right 
of the rightmost root of Cn-i- Consider (jH) at the rightmost root xq of 
Cn-i- As Xq is negative, we know that xq — Xg < 0, and so C„(xo) and 
C4_i(xo) have opposite signs. The claim now follows, since in —00, the 
polynomials Cn{x) and C^_^(x) must converge to the same (infinite) limit 
as their degrees are of the same parity. As C'^_i{x) has no more roots on 
the right of xq, the polynomial Cn{x) must have one. O 

Note that we have in fact proved that the roots of Cn-i{x) and C„(x) 
are interlacing, so the sequence Ci, C2, • • • is a Sturm sequence. 

As an immediate application of the real zeros property, we can determine 
where peak (or peaks) of the sequence Cn,i,Cn,2-, ■ ■ ■ ,Cn,n is- Our tool in 
doing so is the following theorem of Darroch. 
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Theorem 2 141 Let A{x) = X^^^g '^fc^'^ ^ polynomial that has real roots 
only that satisfies A(l) > 0. Let m be an index so that am. = maxo<i<nOi- 
Let /i = A' {1)/A{1). Then we have 

\fi — m\ < 1. 



In particular, a sequence with the real zeros property can have at most 
two peaks. Note that ^'(1) = X]"=o ^(1) = SiLo*^*' therefore 

A'{1)/A{1) is nothing else but the weighted average of the coefficients Oj, 
with i being the weight of Oj. So in the particular case when A{x) = Cn{x), 
we have 

ECn,i 
^' (2n-l)!! 

i 

_ 2n + 1 

where the last step follows from Corollary [H Indeed, (2n-'i)!! j^^* 
probability that a randomly selected Stirling permutation of length n has 
exactly i descents, so sumii ■ (2n-'i)!! J^^^ expected number of descents 
in such permutations. 

Therefore, by Theorem [21 we obtain the following result. 

Theorem 3 Let i be an index so that Cn,i = max^ Cn^k- Then 

1. i = (2n + l)/3 if {2n + l)/3 is an integer, and 

2. i = [(2n + 1)/3J or i = \{2n + l)/3] if (2n + l)/3 is not an integer. 



3 Stirling Permutations and Normal Distribution 

In this section, we prove that the plateaux (equivalently ascents, equiva- 
lently, descents) of Stirling permutations are normally distributed. Our main 
tool is the following result of Bender. Let X„ be a random variable, and let 
an{k) be a triangular array of non-negative real numbers, n = 1, 2, • • • , and 
1 < k < m{n) so that 
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Set <7n(x) = E^i?^'n(^)^'• 

We need to introduce some notation for transforms of the random vari- 
able Z. Let Z = Z - E{Z), let Z = Z/^Var(Z), and let Z„ ^ iV(0, 1) 
mean that Z„ converges in distribution to the standard normal variable. 

Theorem 4 ^2] Let Xn and gn{x) he as above. If gn{x) has real roots only, 
and 

an = \J Yar(Xn) oo, 

then Xn A^(0, 1). 

See [3] for related results. 

We want to use Theorem[5]to prove that the plateaux of permutations in 
Qn are normally distributed. Because of Theorem[Tl all we need for that is to 
prove that the variance of the number of these plateaux converges to infinity 
as n goes to infinity. We will accomplish more by proving an explicit formula 
for this variance. In order to state that formula, let Yn^i be the indicator 
random variable of the event that in a, rand-omly sGlcctcd. Glcment of 
the two copies of i are consecutive, that is, they form a plateau. Note that 

P{Yn,n = 1) = E{Yn,n) = 1- Set Yn = EUl ^n,- 

Theorem 5 For all positive integers n, the equality 

In'^ — 9 

Var(Yn) = — (6) 

^ ^ 18n-9 ^ ^ 

holds. 

Proof: We are going to use the identity Var(y„) = E{Y^) - £'(y„)^. We 
have seen in Corollary □ that E{Yn) = Let s„ = E{Y^). The key 

element of our computations is the following lemma. 

Lemma 1 For all positive integers n, the equality 

2n - 1 4n + 4 

Sn+l = - — — • Sn + — - — . 7 
2n + 1 6 

holds. 

Proof: In order to prove ([7D, we need the following simple facts. 
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Proposition 2 1. For all positive integers n, and all indices i ^ j that 
satisfy 1 < i-, j < n, the equality 



2n- 1 
2n+ 1 

holds. 



2. For all positive integers n and all indices 1 < i < n, the equality 

2n 

holds. 

3. For all indices i < n + 1, the equality 

holds. In particular, E(Yn+i^n+i) = 1- 
Proof: 

1. In order to get an element of Qn+i in which i and j are both plateaux, 
take an element of Qn in which i and j are both plateaux, and insert 
two consecutive copies of n + 1 into any of the 2n — 1 available places, 
that is, anywhere but between the two copies of i or the two copies of 
j- 

2. In order to get an element of Qn+i in which i is a plateau, insert two 
consecutive copies of n + 1 into any of the 2n available slots, that is, 
anywhere but between the two copies of i. 

3. Obvious since n + 1 is always a plateau in elements of Qn+i- 

O 

We return to proving Lemma [H 

Note that s„+i = Z]i<ij<n+i ^he latter can be split 
into partial sums based on whether i or j are equal to n + 1 as follows. 

Sra+l = ^ E{Yn+l^n+lYn+lj) + ^ -E'(^n+l,i^n+l,n+l) 
l<j<«+l l<i<n 

+ Yl E{Yn+l,^Yn+lJ). 
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Based on part 3 of Proposition [21 this simplifies to 

i<j<n+l l<i<n l<i,j<n 

+ Y E{Yn+l,i). 
l<i<n 

Now note that the first sum on the right-hand side is just E(Yn+i)^ the 
second sum is E{Yn+i — ^n+i,n+i) = E{Yn+i) — 1, use part 1 of Proposition 
[2] on the third sum, and part 2 of Proposition [2] on the fourth sum to get 

2n — 1 2n 

Sn+l = 2E{Yn) - 1 + {Sn " E{Yn)) + ^^^^(^n)- 

Recalhng from Corollary [T] that E{Yn) = this reduces to d?]). O 

Using the recursive formula proved in Lemma [H it is routine to prove 
that 

Sn = E{Y„) = . (8) 

\ n) i8n-9 ^ ' 

Therefore, Var(y„) = s„ - E{Ynf' = as claimed. O 



Theorem 6 The distribution of the number of plateaux of elements of Qn 
converges to a normal distribution as n goes to infinity. That is, Yn — > 
iV(0,l). 

Proof: Let X„ = y„, and let gn{x) = (2ra-i)!! ^"(^)- Then Theorem [T] and 
Theorem [5] show that the conditions of Theorem [5] are satisfied, and the 
claim follows from Theorem 31 O 



4 Remarks 

Corollary [T] shows that E{Yn) = (2n + l)/3. It is not difficult to prove that 
E(Yn,n-i) = n}=i 2ra-2jH-i linearity of expectation this proves the 

interesting identity 



where the empty product (indexed by i = 0) is considered to be 1. 

The proof of the equidistribution of the descent and plateau statistics 
we gave is very simple, but it is of recursive nature. A direct bijective proof 
has recently been given by Hyeong-Kwan Ju [5]. 
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